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Abstract 

Mammary cancer is one of the most common neoplasms of dogs, primarily bitches. While studies have been carried out 
identifying differing risk of mammary neoplasia in different dog breeds, few studies have reported associations between 
dog breeds and clinical features such as number of neoplastic lesions found in an individual case or the likelihood of lesions 
being benign or malignant. Such epidemiological studies are essential as a foundation for exploring potential genetic drivers 
of mammary tumour behaviour. Here, we have examined associations between breed, age and neuter status and the odds of 
a diagnosis of a mammary epithelial-origin neoplastic lesion (as opposed to any other histopathological diagnosis from a 
biopsied lesion) as well as the odds of a bitch presenting with either a single mammary lesion or multiple lesions, and the 
odds that those lesions are benign or malignant. The study population consisted of 129,258 samples from bitches, including 
13,401 mammary epithelial neoplasms, submitted for histological assessment to a single histopathology laboratory between 
2008 and 2021. 

In multivariable analysis, breed, age and neuter status were all significantly associated with the odds of a diagnosis of a 
mammary epithelial-origin neoplastic lesion. Smaller breeds were more likely to receive such a diagnosis. In cases diagnosed 
with a mammary epithelial neoplasm, these three factors were also significantly associated with the odds of diagnosis with a 
malignant lesion and of diagnosis with multiple lesions. Notably, while neutered animals were less likely to have a mammary 
epithelial neoplasm diagnosed, and were less likely to have multiple neoplasms, they were more likely to have malignant dis- 
ease. Exploration of the patterns of risk of developing malignant disease, or multiple lesions, across individual breeds showed 
no breed with increased odds of both outcomes. Breeds with altered odds compared to the Crossbreed baseline were either at 
increased risk of malignant disease and decreased risk of multiple lesions, or vice versa, or they were at significantly altered 
odds of one outcome with no change in the other outcome. Our analysis supports the hypothesis that age, neuter status and 
intrinsic biological and genetic factors all combine to influence the biological heterogeneity of canine mammary neoplasia. 
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Introduction 


The mammary gland is a common site for neoplasia in 
dogs [1]. For example, an insurance register study of 
female dogs in Sweden (where only a small proportion 
of bitches, approximately 7%, are neutered) [2] reported 
an overall incidence of 111 cases (including both benign 
and malignant neoplasms) per 10,000 dog years at risk 
(DYAR). Incidence increased with age and varied by breed 
[3]. A study of cancer rates in dogs insured with a UK 
pet insurance company (in the UK approximately 60% of 
female dogs are neutered) [4] reported an estimated age- 
standardised incidence of 205 cases (benign and malig- 
nant) / 100,000 dogs / year [5]. The latter study did not 
distinguish between male and female dogs. 

As with human breast cancer, mammary cancer is more 
common in female than in male dogs [6, 7] and exposure 
to female hormones is likely important for the aetiology of 
mammary neoplasia in bitches, with neutering suggested 
to reduce risk of occurrence [8-11] (D. Varney, D. O'Neill, 
M. O'Neill, D. Church, A. Stell, S. Beck, M. Smalley and 
D. Brodbelt: The epidemiology of mammary tumours in 
bitches under veterinary care in the UK in 2016, submit- 
ted). A recent systematic review questioned the strength 
of the epidemiological evidence supporting an association 
between neutering and mammary tumour risk in dogs, in 
particular relating to the timing of neutering [12], how- 
ever, the data are overall consistent that neutering (particu- 
larly at a younger age) reduces mammary tumour risk [13]. 
If an entire animal develops a mammary tumour, neutering 
combined with mastectomy significantly increases survival 
compared to neutering alone [14] although estrogen may 
have a protective effect in animals with estrogen receptor 
(ER) negative tumours [15], so neutering may be coun- 
terproductive in these cases. Like human breast cancer, 
canine mammary tumours which lack hormone receptors 
tend to be more aggressive [16-18]. Other features of 
similarity between canine and human mammary tumours 
include their age of onset (usually middle-aged to elderly 
individuals), shared environmental exposures, the oppor- 
tunity and common motivation to treat with curative intent 
and their potential for metastasis via the lymphatic system. 
Furthermore, the presence of histological subtypes with 
prognostic significance [19-22] make canine mammary 
tumours an ideal model system within which to investigate 
mammary tumour heterogeneity. 

The key strength of the dog as a disease model, how- 
ever, lies primarily in the genetics of the current breed 
structure of dogs. Modern dog breeds have been highly 
selectively inbred to accentuate particular and desired phe- 
notypic and/or behavioural features. In the process, many 
individual breeds have become predisposed to particular 
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diseases or pathological conditions, either as a direct result 
of the phenotypic feature being selected for (e.g. brachy- 
cephaly and skin fold dermatitis) [23] or because a risk 
allele(s) for the pathology is/are in linkage disequilibrium 
with the genetic variants responsible for the selected fea- 
tures (e.g. Dachshunds and Lafora Disease). Differential 
breed-associated disease predisposition also affects cancer 
risk in dogs [2, 5, 7, 9, 24-29]. For example, we recently 
demonstrated that large dogs are at significantly higher 
risk of osteosarcoma, while brachycephalic dogs (with flat 
faces, such as bulldogs) are at significantly lower risk [30]. 
Breed associations with mammary cancer risk have also 
been demonstrated [3, 7, 9, 24-27]. We have also recently 
identified an increased risk of mammary neoplasia in Eng- 
lish Springer Spaniels, Lhasa Apso and Staffordshire Bull 
Terriers (D. Varney, D. O'Neill, M. O'Neill, D. Church, A. 
Stell, S. Beck, M. Smalley and D. Brodbelt: The epidemi- 
ology of mammary tumours in bitches under veterinary 
care in the UK in 2016, submitted). The increased risk of 
mammary cancer in English Springer Spaniels has previ- 
ously been noted [3, 25], with a suggested link to BRCA/ 
and BRCA2 polymorphisms [31]. 

However, there is little information on whether the breed 
of dog in which a mammary neoplasm develops influences 
the biology of that tumour, over and above any effect of the 
risk of that tumour forming in the first place. If such breed- 
biology associations were elucidated, they would validate 
canine mammary tumours as a model system for investigat- 
ing the genetic drivers of mammary tumour heterogeneity. 

Dogs with mammary neoplasia can present with single 
palpable masses or multiple masses (although anecdotally 
‘single masses’ can be found to be in fact composed of 
more than one lesion upon histological examination). The 
number of lesions found, and their classification as benign 
or malignant, can determine how aggressively the animal 
is clinically managed. Surgical interventions may range 
from simple lumpectomy of a benign adenoma to complete 
removal of the mammary chain, sometimes bilaterally and 
in conjunction with ovario-hysterectomy [1]. Here, we have 
investigated associations between epidemiological risk fac- 
tors (breed, age and neuter status) with two pathological fea- 
tures of canine mammary tumours submitted for histology: 
(1) the number (single or multiple) of lesions on the submit- 
ted mammary chain and (ii) the histotypes (benign or malig- 
nant) reflected by those single or multiple lesions. Data were 
generated from a large archival record of 129,258 biopsy 
samples with associated clinical data from female domestic 
dogs. Samples were submitted to a single diagnostic centre 
between 2008 and 2021, and 13,401 had a confirmed diag- 
nosis of at least one epithelial-origin mammary neoplasm. 

We find a significant association in the current dataset, 
between bitch breed, age and neuter status and the odds 
of being diagnosed with a mammary epithelial neoplasm 
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as opposed to a diagnosis for any other condition. Further- 
more, breed, age and neuter status were also significantly 
associated with whether a case presented with only a sin- 
gle mammary epithelial neoplastic lesion or with multiple 
lesions, as well as with whether benign or malignant disease 
was diagnosed. Of particular note, neutered animals had a 
small but statistically significant increase in odds of diag- 
nosis with a malignant mammary epithelial neoplasm. Our 
findings will form the basis to develop testable hypotheses 
addressing the genetic drivers of mammary tumour biologi- 
cal heterogeneity. 


Methods 
Study Samples and Data Processing 


The study was approved by Cardiff School of Biosciences 
Research Ethics Committee (20-05-01). Information on 
samples submitted to VPG Histopathology, Bristol, UK 
(https://vet.synlab.co.uk/histopathology/) between 2008 
and 2021, with consent for research use, were retrieved by 
a search of the VPG Histopathology block archive database 
for all records listed as canine in origin. In total, 300,592 
records were downloaded as a Microsoft Excel spreadsheet 
with the following information: case number, breed, age at 
sample submission (years and months), sex and neuter sta- 
tus, histological diagnosis. The latter was a semi-standard- 
ised free-text field which contained in some cases a single 
diagnosis, with information on severity and location, whilst 
in other cases there were multiple diagnoses of samples from 
the same case, which might have come from the same or dif- 
ferent anatomical locations. Therefore, each sample record 
was a record of an individual submission from a single ani- 
mal, but might have included biopsies from multiple loca- 
tions. We cannot exclude that, given the extended period of 
time over which samples were submitted, an animal might 
have had samples submitted for analysis more than once. 
Furthermore, as different pathologists had reported on cases 
over the 11-year period, terminology was not uniform even 
for identical lesions. 

Data were cleaned in Microsoft Excel (2013, Microsoft 
Corp.), and Rstudio™ using the packages listed in the R 
script (available at https://github.com/ge8793/Mammary_ 
Public_Data/blob/main/Mammary_Public_Code_Nov22.R). 
Records with duplicate case numbers were identified and, 
where the duplicate was clearly the same record entered 
twice, the entry with the least information (typically lacking 
diagnostic information) was deleted. Where duplicates were 
clearly different submissions which had been given the same 
number in error, both were kept but one was given a letter as 
a prefix to the sample number to make it uniquely identifia- 
ble. Any submission for which either the breed or diagnostic 
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information was absent or ambiguous was excluded. Any 
submissions from wild canids were excluded. Any submis- 
sions with the histological information not in English were 
excluded. Finally, any submission from a male dog or a sub- 
mission for which the sex was unknown was excluded. That 
left 129,258 submissions from female dogs where breed 
information and histology/location of their lesion(s) was 
known. See Table S1 for the details of all submissions. 

Breed descriptions were cleaned to remove typographical 
errors and standardise terminology as detailed previously 
[30]. Some similar breeds were grouped together under 
‘unspecified/other’ both to simplify analysis and allow for 
the fact that full breed designations had not always been 
recorded (e.g. ‘Griffon Bruxellois’ and ‘Basset Griffon 
Vendeen’ were grouped as ‘Griffon (unspecified/other)’; 
‘Springer Spaniel’ and ‘English Springer Spaniel’ were 
grouped as ‘Springer Spaniel (unspecified/other)’). Finally, 
breeds with < 130 individuals (i.e. 0.001%) of the total, were 
grouped together as ‘Other Purebred’. All ‘designer’ crosses 
(purposely bred crossbreeds with contrived names generated 
from two or more purebred breed terms) were included in 
the Crossbreed group. The cleaned and standardised breed 
data comprised the breed variable. 

Neuter status was recorded as female entire, female 
neutered or female, neuter status unknown. Age was 
taken as the date of sample submission and then was cat- 
egorised as: <3.0, 3.0 to<6.0, 6.0 to<9.0, 9.0 to< 12.0 
and > 12.0 years. Bitches for which there was no age data 
recorded were maintained in the analysis, meaning that con- 
fidence intervals and p-values were also generated for the 
‘age unrecorded’ category.. 

Following the initial clean-up, the diagnosis informa- 
tion was searched for the term ‘mammary’. This identified 
13,708 cases in which mammary tissue had been submitted 
for diagnosis (in some cases alongside samples from other 
tissues; Table S2). 

Next, we developed a Natural Language Processing script 
in Python using the Jupyter Notebooks environment on the 
Kaggle platform (https://www.kaggle.com/) to extract diag- 
nostic terms from the records of these 13,708 samples in an 
automated, standardised manner and export these to a CSV 
file. This script was based on the conditional random fields 
approach [32, 33] and also MetaMap (https://Ihncbc.nlm. 
nih.gov/ii/tools/MetaMap.html) to extract key characteristics 
and map these to the National Library of Medicine's Unified 
Medical Language System [34] Metathesaurus. The script is 
available as a Supplementary Data File. 

The extracted histological diagnoses underwent addi- 
tional clean-up using Excel tools to standardise histological 
terminology according to the Surgical Pathology of Tumours 
of Domestic Animals [35] system, first into ductal and non- 
ductal tumours and then non-ductal tumours were further 
classified as simple, non-simple and special types (Table $3) 
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as well as intermediate forms (e.g. carcinoma, tubulopapil- 
lary to solid). Extracted data were also reviewed to ensure 
that where a case had been diagnosed with multiple mam- 
mary lesions of epithelial origin, all had been captured. At 
this stage, a further 307 cases were excluded (Table $1: ‘Tri- 
aged out’) for one of several reasons: no neoplastic lesion 
had been identified (e.g. only mammary hyperplasia); only 
carcinoma in situ was reported (the status of carcinoma in 
situ in the dog is unclear but current guidelines suggest 
it should not be classified as a tumour) [35]; a mammary 
lymph node with metastatic carcinoma was reported but the 
primary site was unknown; the only neoplasm in the mam- 
mary tissue was not of epithelial origin (e.g. fibrosarcoma; 
note, where a sarcoma was recorded as arising in a benign 
mixed tumour, the case was included as a benign epithelial 
lesion and not triaged). This left 13,401 cases with a stand- 
ardised diagnosis of one or more mammary epithelial neo- 
plasms and 115,550 samples with a diagnosis (of any con- 
dition) involving a location other than the mammary gland. 
Table S4 gives the categorisation of all cases as mammary / 
not mammary / triaged out. 

Mammary epithelial cases were then categorised as hav- 
ing malignant lesions (which could include only malignant 
lesions or both benign and malignant lesions) or benign 
lesions and as having only a single diagnosed lesion or 
multiple lesions (which could include multiple mammary 
epithelial neoplastic lesions of the same type or different 
types of mammary epithelial neoplastic lesion). Table $5 
gives the categorisation of the mammary cases as benign/ 
malignant and single/multiple as well as the breakdown of 
the histology of each case. 

Finally, data validation back to the original downloaded 
dataset was carried out to ensure that the diagnostic informa- 
tion and categorisation was linked to the correct submission 
number and signalment data. 


Statistical Analysis 


It was not possible to carry out an accurate power calcula- 
tion a priori for this study, because there were no previous 
studies with which to estimate the effect sizes associated 
with breed and the likelihood of single versus multiple 
mammary epithelial tumours in dogs. Therefore, we chose 
to retrospectively determine the power of the study for one 
breed, Cocker Spaniel (Unspecified/Other) (1001 cases diag- 
nosed with mammary epithelial neoplastic lesions), which 
we found to have a small but still significant increased odds 
of developing multiple lesions compared to Crossbreeds in 
multivariable analysis (OR 1.37, 95%CI 1.24 — 1.84) and one 
breed, Setter (Unspecified/Other) with very few cases (21) 
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but which had the highest odds in this analysis (OR 2.55, 
95%CI 1.06 — 6.11). Using the tool available at https://www. 
openepi.com/Power/PowerCC.htm [36] with methodology 
from Edmunds and colleagues [30, 37], we calculated that 
for Cocker Spaniel (Unspecified/Other) we had 98% power 
while for Setter (Unspecified/Other) we had 50—60% power. 
Obviously, where breeds have fewer representatives in the 
dataset, or effect sizes are smaller, the likelihood of a type-I 
error is increased. 

Statistical methodology was derived from previous stud- 
ies which quantified associations between demographic risk 
factors and cancer in dogs [30, 38] (Varney D, O'Neill D, 
O'Neill M, Church D, Stell A, Beck S, Smalley M, Brodbelt 
D: The epidemiology of mammary tumours in bitches under 
veterinary care in the UK in 2016, submitted). Associations 
between risk factors (breed, age, neuter status) and histolog- 
ical diagnosis of a mammary lesion, and, of the mammary 
tumours only, associations with the development of multi- 
ple mammary epithelial neoplastic lesions or a diagnosis 
of malignant mammary epithelial neoplastic disease were 
determined using binary logistic regression modelling (glm- 
logit function, R-stats package) [30, 39]. A global p-value 
for each was calculated using ANOVA to compare the uni- 
variable regression model with a null model derived from 
the proportions of mammary neoplastic diagnoses, multiple 
mammary epithelial neoplastic lesions and malignant mam- 
mary epithelial neoplastic disease in a comparator group 
(Crossbreeds for the breed variable; age <3 years for the 
age variable; entire animals for the neuter status variable). 
The Imtest package was used to generate global p-values 
in multivariable modelling [30, 40]. Risk factors that were 
significantly associated with the outcomes in univariable 
modelling (p <0.20) were taken forward in multivariable 
modelling that used a manual backwards stepwise elimina- 
tion approach [30, 41]. Final model statistical significance 
was set at the 5% level. The quality of the final multivari- 
able model for each outcome was determined using the area 
under the ROC curve [30, 41]. Pairwise interactions were 
not evaluated for all variables in the final models but instead 
evaluation for interaction was restricted to variables deemed 
to have a relevant biological interaction. Descriptive statis- 
tics from the analyses are given in Table S6 and all results 
are provided in Table S7. 

The association between odds of a mammary neoplasia 
being reported in a breed and the mean body weight for 
a female of that breed (taken from VetCompass estimates 
of average breed mass) [30, 42] (https://www.rve.ac.uk/ 
VetCOMPASS) was explored with simple linear regression 
in GraphPad Prism (GraphPad Software LLC). Analysis 
of distribution of non-simple versus simple tumours was 
carried out using Fisher’s exact test in GraphPad Prism on 
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the distribution of the independent categorical variables 
between single and multiple mammary epithelial neoplastic 
lesion cases. 


Results 
Features of the Study Analysis Sample 


The study population consisted of 129,258 female dogs 
with a sample submitted for routine histological diagnosis 
of any lesion between 2008 and 2021. Of these, 13,401 
cases had received a diagnosis of one or more mammary 
epithelial neoplasms (in some animals, with other addi- 
tional pathology not discussed here) and 115,550 cases 
with a diagnosis (of any condition including both neoplas- 
tic and non-neoplastic disease) involving a location other 
than the mammary gland. Among the study population, the 
most common breeds overall were Crossbreed (22,660 sub- 
missions; 17.5%), Labrador Retriever (16,538 submissions; 
12.8%), Staffordshire Bull Terrier (6,443 submissions; 
5.0%), Cocker Spaniel (6,401 submissions; 4.9%), Jack 
Russell Terrier (5,879 submissions; 4.5%) and Springer 
Spaniel (5,581 submissions; 4.3%) and the median age of 
sample submission was 8 years and 2 months (interquar- 
tile range 6 years — 10 years 3 months). 5.2% of animals 
were entire, 61.4% were neutered, while neuter status was 
not recorded for 33.4% of animals (Tables $1, $4, S6 and 
Fig. 1A). 

The most common breeds diagnosed with a mammary 
epithelial neoplasm were Crossbreed (1,979 cases; 14.8%), 
Labrador Retriever (1,288 cases; 9.6%), Jack Russell Terrier 
(1,123 cases; 8.4%), Cocker Spaniel (1,001 cases; 7.5%), 
Springer Spaniel (684 cases; 5.1%) and Yorkshire Terrier 
(531 cases; 4.0%). The median age of sample submission 
was 9 years (interquartile range 7 years 4 months — 10 years 
9 months). 12.5% of animals were entire, 38.5% were neu- 
tered, while neuter status was not recorded for 49.0% of ani- 
mals (Tables $2, S5, S6 and Fig. 1A). 

The most common breeds diagnosed with a non- 
mammary condition were Crossbreed (20,681 submissions; 
17.9%), Labrador Retriever (15,250 submissions; 13.2%), 
Staffordshire Bull Terrier (5,961 submissions; 5.2%), 
Cocker Spaniel (5,400 submissions; 4.7%), Springer 
Spaniel (4,897 submissions; 4.2%) and Jack Russell 
Terrier (4,756 submissions; 4.1%). The median age 
of sample submission was 8 years (interquartile range 
5 years 8 months — 10 years 3 months). 4.4% of animals 
were entire, 64.1% were neutered, while neuter status was 
not recorded for 31.6% of animals (Tables S1, $4, S6 and 
Fig. 1A). 
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Risk Factors Associated with Diagnosis 
of aMammary Epithelial Neoplasm 


In univariable analysis, there was a significant association 
of breed (p < 0.0001), age (p <0.0001) and neuter status 
(p < 0.0001) with a diagnosis of a mammary epithelial neo- 
plasm versus any other disease, following submission of a 
sample for analysis (Table S7). All factors were therefore 
taken forward into a multivariable model and again all fac- 
tors were highly significant (all p< 0.0001) (Table S7). In the 
multivariable analysis, older age was significantly associated 
with higher odds of a mammary epithelial neoplasm diagno- 
sis (3 to<6 years old: Odds Ratio, OR 13.55, 95%CI 10.47 
— 17.54; 6 to<9 years: OR 32.09, 95%CI 24.92 — 41.32; 9 
to<12 years: OR 36.70, 95%CI 28.51 — 47.26; > 12 years: 
OR 29.33, 95%CI, 22.69 — 37.91). Animals which were neu- 
tered or for which neuter status was unspecified had signifi- 
cantly lower odds of a mammary epithelial neoplasm diag- 
nosis compared to entire bitches (OR neutered 0.17, 95%CI 
0.16 — 0.18; OR neuter status unspecified 0.47, 95% CI 0.44 
— 0.5) (Table S7 and Fig. 1B). 

Breeds including Miniature Dachshund (OR 4.77, 95%CI 
3.74 — 6.09), Pomeranian (OR 3.2, 95%CI 2.39 — 4.29), 
Yorkshire Terrier (OR 3.12, 95%CI 2.78 — 3.49), Toy Poodle 
(OR 2.95, 95%CI 2.13 — 4.09) and English Setter (OR 2.82, 
95%CI 2.08 — 3.83) had significantly higher odds of having 
a diagnosis of a mammary epithelial neoplasm than another 
disease, while Rough Collie (OR 0.33, 95%CI 0.15—0.71, 
Flat-Coated Retriever (OR 0.29, 95%CI 0.19—0.43), St. 
Bernard (OR 0.19, 95%CI 0.06—0.61), Pug (OR 0.18, 
95%CI 0.12—0.29) and Mountain Dog (OR 0.1, 95%CI 
0.04—0.25) had significantly lower odds of having such a 
diagnosis compared to Crossbreeds (Table S7 and Fig. 1B). 

We noted that many of the higher risk breeds were smaller 
dogs while many of the lower risk breeds were larger dogs, 
although there were obvious exceptions (such as the Eng- 
lish Setter and Pug). Individual level body weight data was 
not available in the current study. We therefore plotted the 
breed average body weight for females of each breed (see 
Methods) (Table S8) against the odds ratio of a mammary 
epithelial neoplasm diagnosis for that breed. There was a 
significant inverse correlation (p < 0.001) between the body 
weight and the odds of a diagnosis of mammary epithelial 
neoplasia (Fig. 2A). 


Features of the Cases Diagnosed with a Mammary 
Epithelial Neoplasm 


Of the 13,401 female dogs in our sample diagnosed with one 


or more mammary epithelial neoplasms, 4544 cases (33.9%) 
were diagnosed with a single benign neoplastic epithelial 
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Fig. 1 Presentation of canine mammary neoplasia. A Age distribution 
of all cases in study population (top), non-mammary cases (middle) 
and cases diagnosed with mammary epithelial neoplasia (bottom). B 
Forest plot showing results of multivariable analysis of breed, neuter 
status and age relationship of study population to odds of a diagnosis 
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B Relative Risk of Diagnosis 
with a Mammary Epithelial Neoplasm 
(multivariable analysis) 
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of a mammary epithelial neoplasm versus any other histological diag- 
nosis of disease. Odds ratios and 95% confidence intervals are shown. 
Variables showing a significant difference from base reference value 
(Crossbreed for breed, Entire for neuter status,<3 years for age) are 
indicated in red 
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mammary lesion. 2,314 cases (17.3%) were diagnosed with 
multiple benign lesions. 4,485 cases (33.5%) were diagnosed 
with a single malignant lesion while 2,058 cases (15.3%) 
were diagnosed with multiple lesions including one or more 
malignant neoplasms (Table S7 and Fig. 2B). 

The most commonly diagnosed pathologies were benign 
mixed tumour (4,237 diagnoses), adenoma, complex (3,303 
diagnoses), carcinoma, complex (1,714 diagnoses), carci- 
noma, simple (1,531 diagnoses), carcinoma, ductal (1,139 
diagnoses) and adenoma, simple (928 diagnoses) (Table 1 
shows the complete list of diagnoses). Note that in cases 
presenting with multiple lesions, each case could have more 
than one of these diagnoses and the exact number of lesions 
of each individual type was not necessarily counted (e.g. 
‘multiple simple adenomas and simple carcinomas’). As a 
result, the overall proportions of different diagnoses can- 
not be calculated across the whole sample set. However, 
for cases in which only a single lesion was reported, the 
proportions of the different tumour types can be determined. 
In such single lesion cases, 25.55% (2,226) of tumours were 
benign mixed tumours, 14.84% (1,340) were complex ade- 
nomas, 11.39% (1,028) were simple carcinomas, 10.74% 
(740) were complex carcinomas, 8.00% (722) were ductal 
carcinomas, with all other histotypes being less than 5% of 
the total. The full list is given in Table 2. 


Risk Factors Associated with Clinical Features 
of Mammary Tumours 


In univariable analysis, there was a significant association 
of breed (p <0.0001), age (p < 0.0001) and neuter status 
(p <0.0001) with odds of mammary tumour classification 
as malignant (Table S7). All three were therefore taken for- 
ward into a multivariable model, where again all three were 
significantly associated with risk of malignant disease. 

For age (p < 0.0001 overall), dogs presenting at older ages 
had significantly higher odds of a diagnosis of malignant 
disease relative to younger dogs (6 to<9 years: OR 2.06, 
95%CI 1.16—3.65; 9 to< 12 years: OR 3.12, 95% CI 1.76— 
5.52; > 12 years: OR 4.77, 95%CI 2.68—8.48). For neuter 
status (p<0.0001 overall), neutered animals were signifi- 
cantly more likely to have a diagnosis of malignant disease 
relative to entire bitches (OR 1.25, 95%CI 1.11—1.4) (Table 
S7 and Fig. 3). 

Several breeds had significantly higher or lower odds of a 
diagnosis of malignant mammary epithelial disease relative 
to Crossbreeds in the multivariable analysis. In particular, 
higher odds (ORs from 4.69 to 1.38, p<0.05) were seen in 
the Alaskan Malamute, Samoyed, Newfoundland, Husky, 
Weimaraner, Basset Hound, Cavalier King Charles Spaniel, 
Doberman, Greyhound, German Shepherd Dog, West High- 
land White Terrier, Beagle, Golden Retriever, Bichon and 
the Staffordshire Bull Terrier. Lower odds (ORs from 0.78 to 
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0.32, p<0.05) were seen in the Jack Russell Terrier, Cocker 
Spaniel, Shih Tzu, Chihuahua, Miniature Dachshund, York- 
shire Terrier, Toy Poodle, Irish Setter, Setter and the Tibetan 
Terrier (Table S7 and Fig. 3). 

There was a significant association in univariable analy- 
sis of breed (p <0.0001), age (p < 0.0001) and neuter status 
(p < 0.0001) with a case being diagnosed with multiple ver- 
sus single mammary epithelial neoplastic lesions (Table S7). 
All three were therefore taken forward into a multivariable 
model. All three were significantly associated with altered 
odds of diagnosis of multiple lesions in the multivariable 
model. 

For age (p < 0.0001 overall), dogs presenting at older ages 
had a significantly higher odds of a being diagnosed with 
multiple lesions compared to younger dogs (9 to < 12 years: 
OR 2.42, 95%CI 1.28—4.59;> 12 years: OR 2.23, 95%CI 
1.17—4.25). For neuter status (p < 0.0001 overall), neutered 
animals had significantly lower odds of being diagnosed 
with multiple lesions than entire bitches (OR 0.71, 95%CI 
0.63—0.8) (Table S7 and Fig. 4A). 

A number of breeds had significantly higher or lower 
odds of being diagnosed with multiple lesions relative to 
Crossbreeds in the multivariable analysis. Higher odds (ORs 
from 2.55 to 1.37, p< 0.05) were seen in the Setter, Ameri- 
can Bulldog, Spaniel, Akita, Yorkshire Terrier and Cocker 
Spaniel. Lower odds (ORs from 0.75 to 0.12, p< 0.05) were 
seen in the Border Collie, Staffordshire Bull Terrier, Collie, 
Golden Retriever, Husky, Cavalier King Charles Spaniel, 
Shepherd Dog and Pug (Table S7 and Fig. 4A). 

The changes in ORs of developing malignant / multiple 
neoplasms relative to baseline are summarised for the breed, 
age and neuter status risk factors in Table S9. 


Non-simple Histotypes are Significantly more Likely 
to be Found in Tumours from Cases Presenting 
with Multiple Lesions 


Mammary epithelial neoplastic lesions are classified into 
simple, non-simple, ductal and special histotypes [35]. 
Non-simple neoplasms contain proliferating neoplastic 
populations of more than one lineage, including proliferat- 
ing luminal epithelial and myoepithelial cells and poten- 
tially mesenchymal elements. We speculated whether cases 
which present with multiple neoplasms are more likely to 
have a higher proportion of non-simple neoplasms, com- 
posed of multiple cellular lineages, because of a general- 
ised increased risk of carcinogenesis across all cells in the 
entire mammary chain. Therefore, we assessed whether such 
non-simple lesions were more frequent in animals diagnosed 
with multiple lesions. Considering only cases in which either 
simple or non-simple lesions were diagnosed, either alone 
or in combination with other types (excluding those cases 
in which only ductal lesions, special types or intermediate 
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<Fig.2 A Simple linear regression of mean breed body weight for 
purebred female dogs>18 months old against odds of a diagnosis 
of mammary epithelial neoplasm relative to Crossbreed dogs. Each 
point represents one breed; points in blue/red represent breeds with 
significantly lower (blue) or higher (red) odds. There is a highly sig- 
nificant (py =0.0004) inverse correlation between weight and diagno- 
sis of a mammary epithelial neoplasm. See Table S8 for the details 
for each breed. B Cartoon indicating categories (single benign, mul- 
tiple benign, single malignant, multiple malignant) into which dogs 
presenting with mammary tumours can be divided. Blue circles rep- 
resent benign tumours, red cells malignant tumours. White ellipses 
indicate the five pairs of mammary fat pads 


forms were diagnosed) (11,459 cases), non-simple lesions 
were identified significantly more frequently (89.74% of 
cases) when multiple lesions were diagnosed compared 
to cases with only single lesions (66.48% of which were 
non-simple) (Fisher’s exact test of categorical variables, 
p <0.0001) (Fig. 4B). 


Discussion 


Here, we have analysed histological diagnoses from an 
archive of 129,258 records of submissions of samples from 
bitches, including 13,401 mammary epithelial neoplasms, 
sent for histological assessment to a single histopathology 
laboratory between 2008 and 2021. In multivariable analy- 
sis, we found a significant association between breed, age 
and neuter status and the odds that a sample submitted for 
histological diagnosis would be diagnosed as a mammary 
epithelial tumour, as opposed to any other diagnosis. We 
also found a significant association in multivariable analysis 
between breed, age and neuter status and odds of diagnosis 
of a malignant versus benign mammary epithelial neoplasm 
and also the odds of being diagnosed with multiple versus 
single mammary epithelial neoplasms. Our study was not 
designed to determine how breed, age and neuter status 
affect odds of mammary tumour diagnosis relative to the 
general canine population as this requires comparison to 
the entire population of canines in a particular region as 
a denominator. We (D. Varney, D. O'Neill, M. O'Neill, D. 
Church, A. Stell, S. Beck, M. Smalley and D. Brodbelt: The 
epidemiology of mammary tumours in bitches under veteri- 
nary care in the UK in 2016, submitted) and others [3, 7, 9, 
24-27] have, however, addressed this question elsewhere. 
Our study was primarily designed to address how these vari- 
ables might influence clinical features of tumours once histo- 
logically diagnosed; however, we also took advantage of our 
dataset to understand how they affect the likelihood that any 
sample undergoing histological assessment was diagnosed 
as a mammary epithelial neoplasm rather than any other 
histological diagnosis. 

An increased likelihood of a sample submitted being his- 
tologically diagnosed as a mammary epithelial tumour in 
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older animals reflects mammary neoplasia (and in particu- 
lar adenoma/ adenocarcinoma) as a disease of older dogs 
[29, 43]. Similarly, an increased likelihood in entire ani- 
mals would be consistent with studies showing that neutered 
animals are protected against mammary neoplasia [1 1-13, 
44-47] (Varney D, O'Neill D, O'Neill M, Church D, Stell A, 
Beck S, Smalley M, Brodbelt D: The epidemiology of mam- 
mary tumours in bitches under veterinary care in the UK in 
2016, submitted). However, caution should be exercised in 
the interpretation of our results, as our findings specifically 
relate to samples submitted for histological examination and 
sample submission bias is likely to be a factor. This is dis- 
cussed further below. 

Using the large dataset, of 13,401 cases submitted with 
one or more confirmed epithelial-origin mammary neo- 
plastic lesions (sometimes in addition to other pathologies) 
available to us, we next asked whether the age, neuter status 
and breed variables were also associated with the clinical 
behaviour of the mammary neoplasms in these cases. We 
examined the odds of developing a malignant (as opposed to 
benign) lesion and the odds of being diagnosed with multiple 
(as opposed to single) lesions. In all cases, our analysis was 
based on odds compared to a baseline group (age <3 years; 
‘entire’ neuter status and Crossbreed animals) (summarised 
in Table S9). 

A number of studies have assessed the relationship 
between age at diagnosis and the risk of being diagnosed 
with benign or malignant mammary lesions [7, 20, 26, 27, 
48-50] and of being diagnosed with single or multiple mam- 
mary lesions [26] (although the proportion of dogs with 
single or multiple lesions varies widely between reports, 
potentially a result of differing proportions of breeds within 
the study populations) [49, 51]. We find in multivariable 
analysis that older dogs are more likely to be diagnosed with 
malignant disease and with multiple lesions. This is consist- 
ent with previous reports that dogs with malignant tumours 
are more likely to be older than dogs with benign tumours 
[26, 48, 49], that age is an independent prognostic factor cor- 
relating with poor survival [20] and that older animals tend 
to have higher grade ER negative tumours [27]. A number 
of mechanisms may drive this increased risk of malignancy 
in older dogs. For example, it may simply be a stochastic 
process, with more mutations required to generate a malig- 
nant tumour than a benign one. An example of this may 
be carcinoma arising in complex adenoma/benign mixed 
tumour. Cellular aging may also be a factor as a result of 
telomere shortening, which would lead to increased risk of 
genomic abnormalities [52]. Finally, older individuals tend 
to have a dysfunctional immune system [53], with potential 
for reduced immunosurveillance permitting tumour growth 
and metastasis [54]. 

Neuter status is also well established as a risk factor for 
development of mammary epithelial neoplasia. Multiple 
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Histotype Number 
of diagno- 
ses 

Benign mixed tumour 4237 
Adenoma, complex 3303 
Carcinoma, complex 1714 
Carcinoma, simple 1531 
Carcinoma, ductal 1139 
Adenoma, simple 928 
Adenoma, ductal, papillary 878 
Adenoma, ductal 692 
Carcinoma, simple, tubular 463 
Carcinoma, mixed 374 
Carcinoma, ductal, papillary 368 
Carcinoma, simple, solid 304 
Carcinoma, simple, tubulopapillary 269 
Carcinoma arising in a complex adenoma/benign mixed tumour 158 
Carcinoma, simple, anaplastic 135 
Adenoma, simple, tubular 133 
Carcinoma, simple, comedo 108 
Carcinosarcoma 70 
Carcinoma, simple, papillary 64 
Adenoma, simple, tubulopapillary 54 
Carcinoma and malignant myoepithelioma 44 
Carcinoma, simple, tubular to solid 43 
Carcinoma, simple, micropapillary 42 
Carcinoma, adenosquamous 33 
Carcinoma, lipid-rich 32 
Malignant myoepithelioma 21 
Adenoma, simple, papillary 20 
Myoepithelioma 17 
Carcinoma, simple, cribriform 16 
Carcinoma, inflammatory 15 
Carcinoma, spindle cell 13 
Carcinoma, simple, tubulopapillary to solid* 8 
Carcinoma, simple, solid to comedo* 6 
Carcinoma, simple, solid to tubular* 5 
Carcinoma, simple, tubular to solid to comedo* 3 
Carcinoma, simple, tubulopapillary to comedo* 3 
Carcinoma, complex to micropapillary* 2 
Carcinoma, mucinous 2 
Carcinoma, simple, solid to anaplastic* 2 
Carcinoma, simple, solid to micropapillary* 2 
Carcinoma, simple, tubular to anaplastic* 2 
Carcinoma, simple, tubulopapillary to anaplastic* 2 
Adenoma, simple to complex to ductal, papillary* 1 
Adenoma, simple, tubular to complex* 1 
Carcinoma, complex to adenosquamous* 1 
Carcinoma, complex to tubulopapillary* 1 
Carcinoma, simple, cribriform to solid to comedo* 1 
Carcinoma, simple, papillary to anaplastic* 1 
Carcinoma, simple, solid to adenosquamous* 1 
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Table 1 (continued) Histotype 


Carcinoma, simple, solid to comedo to anaplastic* 
Carcinoma, simple, solid to complex* 

Carcinoma, simple, solid to cribriform* 
Carcinoma, simple, solid to cribriform to comedo* 
Carcinoma, simple, solid to spindle cell* 
Carcinoma, simple, solid to tubular to comedo* 


Carcinoma, simple, solid to tubulopapillary to comedo 


Carcinoma, simple, solid, cystic 


Carcinoma, simple, tubular to comedo* 
Carcinoma, simple, tubular to cribriform* 
Carcinoma, simple, tubular to mixed* 
Carcinoma, simple, tubular to mucinous* 


Carcinoma, simple, tubular to spindle cell* 


Number 
of diagno- 
ses 


a 


fe ee 


* Intermediate forms 


studies, including our own, have demonstrated that neutered 
bitches have a reduced risk of neoplastic mammary disease 
compared to entire animals [10, 11, 44-47] (summarised in 
[13]). Here, we have now shown in multivariable analysis 
that neutered animals are also at reduced odds of develop- 
ing multiple lesions; however, if they do develop neoplastic 
mammary disease, they are at increased odds that this will be 
malignant. It is notable that ER negative canine mammary 
tumours (CMT) are more likely to be malignant than ER 
positive CMT [16-18]. Furthermore, a study which exam- 
ined the link between neutering, serum estrogen levels and 
CMT hormone receptor status demonstrated that ER nega- 
tive tumours in entire animals with high serum hormone 
levels had a longer time to metastasis than such tumours in 
neutered animals [15]. A protective effect of estrogen via 
non-receptor mechanisms was suggested. 

Some studies have found no links between breed and 
predisposition to develop malignant as opposed to benign 
lesions [50] and no difference in malignancy between cases 
presenting with single or multiple masses [55]. However, 
others have suggested particular breeds are more likely to 
develop malignant mammary tumours (Samoyed, Dober- 
mann, Schnauzer and Yorkshire Terrier) [56] (although 
notably when cancer of all sites was considered Aupperle- 
Lellbach and colleagues found Yorkshire Terriers among 
the breeds more likely to develop benign tumours [24]). 
There is little, if any, information on the link between 
breed and risk of developing single as opposed to mul- 
tiple lesions. We find a significant association between 
breed and both of these aspects of tumour biology. The 
patterns of odds of developing malignant disease, or multi- 
ple lesions, across individual breeds showed that no breed 
was found to be at increased odds of both outcomes. They 


were either at increased odds of malignant disease and 
decreased odds of multiple lesions (e.g. King Charles 
Cavalier Spaniel, Golden Retriever, Husky, Staffordshire 
Bull Terrier) or vice versa (e.g. Cocker Spaniel, Setter, 
Yorkshire Terrier) or they were at significantly altered 
odds of one outcome with no change in the other outcome. 
The simultaneous appearance of multiple mammary 
epithelial neoplasms in cases presenting with multiple 
lesions suggests a ‘field cancerisation’ model. Field can- 
cerisation was first proposed by Slaughter and colleagues 
in 1953, to describe the development, in oral squamous 
cell carcinoma, of regions of the oral epithelium which are 
clinically apparently normal but in which multiple inde- 
pendent primary squamous cell carcinomas continually 
arise [57]. It is now considered that field cancerisation is 
the product of an underlying preneoplastic stem cell which 
has acquired mutations giving it a competitive advantage 
overing neighbouring cells, allowing clonal progeny of the 
original preneoplastic cell to spread. This then creates a 
large target population requiring fewer mutational events 
for full transformation [57-61]. We suggest that the pres- 
ence of multiple independent tumours, a high proportion 
of which are non-simple tumours composed of multiple 
cell lineages, supports a model that in some dogs the entire 
mammary epithelium is a field of preneoplastic (stem) 
cells. Furthermore, we suggest that there can be genetic 
factors predisposing to the development of a preneoplastic 
mammary field. A number of candidate genes have been 
identified which affect cell competition and have known 
roles in cancer, for example p53 [60], so these would be 
an excellent starting point for future studies. 
Interestingly, Gunnes and colleagues [49] reported that 
in bitches presenting with multiple tumours, the chance 
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Table 2 Percentage of tumour histotypes diagnosed in cases with sin- 
gle lesions 


Histotype Number Percent of cases 
of cases 
Benign mixed tumour 2126 23.55% 
Adenoma, complex 1340 14.84% 
Carcinoma, simple 1028 11.39% 
Carcinoma, complex 970 10.74% 
Carcinoma, ductal 722 8.00% 
Adenoma, ductal, papillary 391 4.33% 
Carcinoma, simple, tubular 313 3.47% 
Adenoma, ductal 309 3.42% 
Adenoma, simple 293 3.25% 
Carcinoma, ductal, papillary 234 2.59% 
Carcinoma, mixed 233 2.58% 
Carcinoma, simple, solid 220 2.44% 
Carcinoma, simple, tubulopapillary 180 1.99% 
Carcinoma, simple, anaplastic 102 1.13% 
Carcinoma arising in a complex 89 0.99% 
adenoma/benign mixed tumour 
Carcinoma, simple, comedo 84 0.93% 
Carcinosarcoma 53 0.59% 
Adenoma, simple, tubular 47 0.52% 
Carcinoma, simple, papillary 46 0.51% 
Carcinoma and malignant myoepithe- 33 0.37% 
lioma 

Carcinoma, lipid-rich 30 0.33% 
Carcinoma, simple, micropapillary 27 0.30% 
Adenoma, simple, tubulopapillary 24 0.27% 
Carcinoma, adenosquamous 21 0.23% 
Carcinoma, inflammatory 12 0.13% 
Carcinoma, spindle cell 9 0.10% 
Malignant myoepithelioma 8 0.09% 
Adenoma, simple, papillary 7 0.08% 
Carcinoma, simple, cribriform 7 0.08% 
Myoepithelioma 5 0.06% 
Carcinoma, mucinous 1 0.01% 
Intermediate forms 65 0.72% 


that two tumours would have the same diagnosis and 
level of malignancy was greater than would be expected 
by chance alone, leading them to conclude the existence 
of a hormone-driven field cancerisation effect. They also 
suggest that there might be breed variations in predisposi- 
tion to malignant (as opposed to benign) tumours but their 
study was not sufficiently powered to identify specific at- 
risk breeds. 

Links between canine mammary tumour histological 
subtype and prognosis are well established [19-22] and the 
prevalence of different histological subtypes within canine 
mammary neoplasia has also been assessed in a number of 
studies, for example [7, 49, 50]. Results differ widely from 
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study to study, likely as a result of differing study popula- 
tions and differing histological interpretations. Salas and 
colleagues report similar frequencies of benign and malig- 
nant tumours in their study, with more epithelial-type than 
mixed neoplasms [7]; Gunnes and colleagues reported 61% 
of examined tumours were benign, 39% were malignant, 
with complex adenoma and complex carcinoma the most 
frequent diagnosis in each category [49]; Ariyarathna and 
colleagues report 56% of examined tumours as malignant 
(simple carcinomas being most common) and 44% benign 
(benign mixed tumours being most common). We find that in 
cases with single lesions being diagnosed, the most common 
lesions were the benign mixed tumour (23.6%), adenoma, 
complex (14.8%), carcinoma, simple (subtype not otherwise 
specified (11.4%), carcinoma, complex (10.7%) and carci- 
noma, ductal (8.0%). 

While a diagnosis of mixed and complex (non-simple) 
mammary neoplasms is common in the dog, no studies 
have yet reported a link between a diagnosis of non-sim- 
ple tumours and development of multiple lesions. The link 
between complex/mixed lesions and presentation with mul- 
tiple lesions (with a potential field cancerisation effect) is 
of particular interest as the origin of complex/mixed lesions 
remains unclear. Complex tumours contain distinct prolif- 
erating luminal epithelial and myoepithelial populations 
[35] Mixed tumours contain, in addition, mesenchymal 
elements [35]. It is possible that complex/mixed tumours 
are polyclonal in origin, with separate transformed luminal 
and myoepithelial cells (in the case of a complex tumour) 
or luminal, myoepithelial and mesenchymal stem cells (in 
the case of mixed tumours) all contributing to a neoplasm 
presumably in response to a highly localised tumour promot- 
ing factor, such as an inflammatory signal. Alternatively, 
these tumours might arise monoclonally from a mammary 
epithelial stem cell capable of undergoing both luminal and 
myoepithelial differentiation as well as metaplastic potential. 
Previous studies have addressed this question by analysis of 
mitochondrial DNA mutations in the epithelial and mesen- 
chymal elements [62] or by analysis of immunohistochemi- 
cal staining patterns and DNA ploidy of the different com- 
ponents of the tumour [63]. These analyses have suggested 
some tumours may be polyclonal while others may be mono- 
clonal, but many of the analyses have not been informative. 
A definitive answer on the aetiology of these tumours awaits 
further study. 

Our analysis supports the hypothesis that neuter status, 
age and intrinsic biological and genetic factors all influence 
the heterogeneity of clinical presentation of canine mam- 
mary neoplasia. The term ‘heterogeneity’ includes both 
“intra-tumour heterogeneity’, the genetic, epigenetic, pheno- 
typic and/or behavioural differences in cells within a tumour 
(including both neoplastic cells and non-neoplastic cells 
such as tumour-associated fibroblasts and macrophages), and 
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Fig.3 Correlation of features of study population with risk of malig- > 


nancy of neoplastic epithelial mammary lesions. Forest plot show- 
ing results of multivariable analysis of breed, neuter status and age 
relationship of study population to odds of a diagnosis of a malig- 
nant mammary epithelial neoplasm. Odds ratios and 95% confidence 
intervals are shown. Variables showing a significant difference from 
base reference value (Crossbreed for breed, Entire for neuter sta- 
tus,<3 years for age) indicated in red. See also Table S7. *Odds 
Ratio with 95% confidence intervals cannot be provided for Mountain 
Dog (Unspecified/Other) as all five cases in this breed were malig- 
nant. Confidence intervals are very wide for breeds such as the St 
Bernard as they are very rare breeds in the data set 


inter-tumour heterogeneity, the classification of tumours as 
different histotypes which may have different clinical behav- 
iours (e.g. benign or malignant disease) and approaches to 
therapy [64]. Tumour heterogeneity may be considered to 
arise from two processes. In the early phase of tumour devel- 
opment, it arises from the interaction between cell of origin 
of the tumour and the initiating genetic lesions occurring in 
that tumour [65-67]. Then, as the tumour progresses, the 
random generation of progressively more mutated clones, 
combined with selective pressures on these variant clones, 
lead to tumour heterogeneity through a Darwinian evolution- 
ary process [68]. We would argue that the former process is 
the main driver of inter-tumour heterogeneity while the latter 
is the main driver of intra-tumour heterogeneity. Of course, 
these are not clear-cut divisions, considering that a change 
in the proportion of different cell types within a tumour (for 
example, cells that express hormone receptors in a mammary 
neoplasm) could lead to a breast cancer being reclassified 
as progressing from ER positive to ER negative disease i.e. 
a change in intra-tumour heterogeneity leads to a change in 
inter-tumour heterogeneity. 

Modelling the interaction between the cell of origin 
and genetic lesion as a determinant of mammary tumour 
heterogeneity in genetically modified mice has confirmed 
the principles that tumour histotype is driven by the inter- 
actions between cell of tumour origin, initiating genetic 
lesion and in some cases developmental history of the 
gland [65-67, 69]. However, this is an artificial system 
in which candidate genes are bred into mice in condi- 
tional knockout/overexpression scenarios and in which 
a limited number of tumour histotypes develop with (for 
the most part) little clinical relevance either to human or 
veterinary medicine. In contrast, canine mammary tumours 
combined with the power of dog genetics offer a system 
in which, rather than choosing the genes of interest and 
working forward to understand what, if any, effect they 
have on tumour phenotype, one can work backwards from 
tumour phenotype to elucidate the underlying genetics. The 
first stage of this is to use an epidemiological approach 
to establish associations between tumour biology and 
breed and develop hypotheses. These can then be tested in 
case-control genomic studies within and between breeds 
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Relative Risk of Diagnosis with 
Malignant Mammary Epithelial Disease 
(multivariable analysis) 

Benign Malignant 


Crossbreed 

Alaskan Malamute 

amoyed 

Newfoundland 

Husky (Unspecified/Other) 
Dalmatian 

Kerry Blue Terrier 
Weimaraner 

Basset Hound 

Cavalier King Charles Spaniel 
Italian Spinone 

Doberman 

Greyhound (Unspecified/Other) 
German She ee Dog 

it. Bernard 

Old English Sheepdog 

West Highland White Terrier 
iniature Schnauzer 

Bull Terrier (Unspecified/Other: 
Beagle (Unspecified/Other, 
Soft-Coated Wheaten Terrier 
Retriever (Unspecified/Other 
Bulldog (Unspecified/Other 
Shar Pei 

Mastiff (Unspecified/Other) 
logue de Bordeaux 
Flat-Coated Retriever 
Leonberger 

Golden Retriever 


Shepherd Dog (iemehotdaer 
Bichon (Unspecified/Other, 


hasa Apso 
Pointer Unspecifediotner 
Akita (Unspecified/Other: 
Rhodesian Ridgeback 
Collie (Unspecified/Other) 
Cairn Terrier 


Schnauzer (Unspecified/Other) 
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[70] to identify loci associated with the particular aspect 
of biology (e.g. presentation with multiple tumours or 
a diagnosis of malignant disease) which can be taken 
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<Fig.4 Correlation of features of study population with risk of pre- 
senting with multiple neoplastic epithelial mammary lesions. A For- 
est plot showing results of multivariable analysis of breed, neuter 
status and age relationship of study population to odds of being diag- 
nosed with multiple mammary epithelial neoplasms. Odds ratios and 
95% confidence intervals are shown. Variables showing a significant 
difference from base reference value (Crossbreed for breed, Entire for 
neuter status,<3 years for age) indicated in red. See also Table $7. 
*Odds Ratio with 95% confidence intervals cannot be provided for 
the St Bernard as all three cases in this breed were single cases. B 
Proportion of cases presenting with single and multiple lesions diag- 
nosed with one or more Non-Simple neoplasms, as opposed to only 
Simple neoplasms. ****p <0.0001 (Fisher’s exact test) 


forward into mechanistic studies. The results we present 
here represent the first stage of this process. 

Our study has limitations and caveats, the principle one 
being selection bias (i.e. the dogs in the current study may 
not be representative of the general caseload of canine 
mammary tumours as a whole) [71]. There are many fac- 
tors which affect clinical decision making when deciding 
whether to submit lesions for histological analysis and these 
could all bias the population of dogs that appear in the cur- 
rent study (a ‘biopsy-only’ dataset). For example, dogs with 
very severe disease or poorly resectable lesions may undergo 
euthanasia or palliative care rather than surgical resection 
and histological analysis. In contrast, an entire bitch which 
presents at a primary care practice with a mammary mass 
may be more likely to have that mass sent for histology than 
a neutered bitch, simply because of the previously postulated 
links between neuter status and mammary tumour risk. Fur- 
thermore, financial limitations regarding the cost of surgery 
and biopsy submission may mean that dogs owned by own- 
ers of lower socio-economic status or uninsured dogs may be 
less likely to appear in the current dataset. Previous studies 
have shown that both breed and neutering status in dogs 
are associated with owner socioeconomic status, therefore 
if owner demographics are biased in the dataset, the distri- 
bution of breeds or neuter status may be altered [4]. The 
socio-economic circumstances or age of owners may also 
have resulted in an underrepresentation of benign lesions, for 
example, if individuals with lower economic means or who 
face challenges accessing veterinary care were less inclined 
to take a dog with a mammary mass for veterinary atten- 
tion unless it shows obvious signs of malignancy, such as 
extremely rapid growth or obvious morbidity. Furthermore, 
in certain breeds a single benign lesion may be more difficult 
for an owner to detect owing to body shape or conforma- 
tion, whereas multiple lesions may be more likely to lead 
to an owner seeking attention. We found an inverse correla- 
tion between average breed body weight and likelihood of 
a histopathological sample being diagnosed as a mammary 
epithelial tumour. However, while a predisposition for mam- 
mary tumour development in small dogs has been reported 
[7, 29] population-based case-control studies have identified 
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both small and large breeds as being at higher odds of devel- 
oping mammary epithelial neoplasia [3, 7, 9, 24—26]. It may 
be that smaller dogs are more likely to be picked up and have 
mammary masses discovered by owners. Another alternative 
is that small dogs live longer than large dogs, and mammary 
neoplasia is a disease of older animals. Overall, selection 
bias could have affected the observed associations between 
breed, age and neuter status and mammary tumour number 
or histotype. 

The study may also have been affected by missing data 
because in veterinary clinical practice, there is not a stand- 
ardised approach to lumpectomy vs mammectomy vs full 
mammary strip removal. While full mammary strips were 
frequently submitted (and indeed such submissions were 
encouraged), we cannot exclude that some practices did not 
submit all tissue or all lesions present in a dog for analysis, 
leading to reporting of cases with more than one lesion as 
single lesion cases. However, it is more likely that in cases 
which present with only one, or a small number of, palpable 
lesions, all will be sent for analysis, whereas in cases pre- 
senting with many such lesions only a representative sam- 
ple might be examined. Such an approach would still result 
in the case being correctly categorised as having multiple 
lesions. It is more likely that cases could be incorrectly cat- 
egorised as having single lesions, if only a single palpable 
lesion is detected and sampled, but small, clinically unde- 
tectable, lesions are already present. 

Once tissue has been submitted, given the extensive 
experience of VPG Histopathology, it is highly unlikely that 
cases would be incorrectly categorised as ‘benign’ rather 
than ‘malignant’, even during histological analysis of large 
pieces of tissue (such as whole mammary strips) owing to 
the meticulous nature of the analysis. However, we cannot 
definitively exclude that occult malignant cell clusters may 
be present in an otherwise benign tissue and, if small, these 
cases could theoretically be miscategorised as benign-only. 

We did not carry out Bonferroni corrections for multiple 
comparisons, which may have caused some Type-! errors. 
However, this is usually considered overly stringent where 
there may be correlation between variables [72]. Further- 
more, although unlikely, we cannot exclude the possibility 
that there may be multiple samples submitted for histologi- 
cal examination from the same animal on different occa- 
sions. To definitively exclude this would require access to 
information that would identify owners and therefore such 
data was withheld from the current study. 

We also did not account for the expected breed lifespan 
when considering risk factors for cancer. However, breed 
and neuter status may affect the years-at-risk of dogs, and 
thus certain breeds with longer lifespans, or neutered dogs 
(shown to live longer on average) may appear to be more at 
risk simply because they experience more years-at risk of 
disease. The inclusion of this complex variable was beyond 


va Springer 


6 Page 16 of 18 


the scope of this study because a reliable published lifespan 
could not be sourced for all of the included breeds, however 
future analyses considering years-at-risk would be valuable 
for validation of the conclusions presented here. 

Importantly, our study represents only a snapshot in the 
clinical pathway of each animal diagnosed. The diagnostic 
records available to us are not linked to clinical outcome, 
so we do not know whether a diagnosis of single or multi- 
ple lesions affects prognosis although it is clear that malig- 
nant disease has a worse prognosis than benign disease 
(the prognostic significance of histological subtypes is also 
established) [19, 20]. Both benign and malignant disease 
was seen in cases diagnosed with both single and multiple 
lesions, and previous findings have also shown no relation- 
ship between the presence of single and multiple masses 
and a malignant diagnosis [55]. It is unclear whether having 
multiple lesions is an independent prognostic factor. To test 
this would require a prospective study in order to ensure that 
potential confounding factors are controlled for (e.g. the lack 
of a uniform approach to treating mammary masses). Testing 
this would be an important follow-up to the current study. 

Our study also only represents a snapshot of the underly- 
ing biology of the disease. We cannot exclude, for example, 
that a bitch diagnosed with a single malignant lesion might 
have, in the future, gone on to develop multiple lesions, or 
that an animal with multiple lesions that has a full mammary 
strip, and all those lesions are diagnosed as benign after 
surgery, might have developed a malignant tumour at a later 
date if no interventions were performed. To draw definitive 
conclusions that a case with a single lesion diagnosis would 
not have later been diagnosed with multiple lesions, or that 
a case with one or more benign lesions would not have later 
been diagnosed with one or more malignant lesions would 
require surgeons to perform only individual lumpectomies 
on dogs presenting with mammary masses, no matter how 
many masses they present with and no matter how many 
times they return to clinic. This is not consistent with the 
best welfare of the animals involved. However, as our study 
sample overall captures a very large population over a period 
of years (canine lifespans), generalisations seem reasonable, 
as the specific sampling time in the course of the disease 
becomes less important for each individual. 

Therefore, although the caveats above must be kept in 
mind, we have found significant associations between the 
breed, age and neuter status of a bitch, and whether a pres- 
entation for mammary neoplasia is likely to be for single 
lesions or multiple lesions, and whether those lesions are 
likely to be benign or malignant. Furthermore, non-simple 
lesions are enriched in cases presented with multiple neo- 
plasms. We therefore suggest that underlying genetic fac- 
tors can affect tumour heterogeneity, by influencing clinical 
behaviour (the development of benign or malignant disease), 
tumour number and cellular composition. Environmental 
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influences such as aging are also likely to play a role. The 
mammary epithelium of breeds at higher risk of presenting 
with multiple neoplastic mammary lesions may be a pre- 
neoplastic field genetically primed for tumour development. 
Case-—control genomic studies, and mechanistic evaluation, 
have the potential to identify, in an unbiased manner, genes 
driving mammary tumour behaviours and thus such studies 
could substantially advance our understanding of the drivers 
of mammary tumour formation and tumour heterogeneity 
and, ultimately, identify new targets for therapy. 
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